Hydrodynamics of driven granular gases 
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Hydrodynamic equations for granular gases driven by the Fokker-Planck operator are derived. 
Transport coefficients appeared in Navier-Stokes order change from the values of a free cooling state 
to those of a steady state. 
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I. INTRODUCTION 

The gas kinetic theory of elastic particles played impor- 
tant roles in history of statistical mechanics 1, 2J. When 
there are inelastic interactions among particles, the be- 
havior of collections of particles is completely different 
from that of elastic particles: There are no equilibrium 
states and any spatially homogeneous states are no longer 
stable. Such a collection of the particles having inelastic 
interactions is called the granular gas whose physical real- 
ization can be observed in rings of planets, small planets, 
suspended particles in fluidized beds, aerosols, and rapid 
granular flows 

A typical example of granular gases can be found in 
aerosols or suspensions etc. in which the buoyancy is bal- 
anced with gravity. 0, 0| In dense suspensions the hydro- 
dynamic interaction among particles are important 
but effects of the air may be regarded as a thermostat 
driven by the Langevin forces in dilute suspensions. Most 
of researches for granular gases are interested in undriven 
systems which are difficult to be achieved in actual ex- 
periences. We believe that more systematic studies for 
driven granular gases are required. Montanero and San- 
tos analyzed statistical properties of the homogeneous 
state and hydrodyanmics of granular gases the white 
noise thermostat and the Gaussian thermostat. 0,13 Car- 
rilio et a/.j3| and Pagnani et aZ.^^l analyzed granular 
gases in driven systems by the Langevin force which in- 
cludes both the white noise with the friction force in 
proportion to the velocity of particles, but they are not 
interested in hydrodynamics of such the granular fluid. 
Since we believe that the model driven by such the ther- 
mostat, we call the Langevin thermostat, can describe 
physical situations of dilute suspensions, we need to clar- 
ify properties of the hydrodynamics of such the system. 
The main purposes of this paper is to investigate the ef- 
fects of Langevin thermostat for transport coefficients 

The organization of this paper is as follows. In the 
next section, we introduce our model and the framework 
of the Chapman-Enskog method for the analysis of gran- 
ular gases. In section III, we calculate the valuables such 
as the granular temperature and the fourth cumulant in 
a homogeneous state. In section IV, we obtain the trans- 
port coefficients such as the viscosity and the heat con- 
ductivity. In section V, we discuss and summarize our 
results. 



II. FRAMEWORK 

We consider rarefied gases of smooth identical particles 
with the mass m, the velocity v and the diameter a. The 
distribution function /(r,v, i) in our system obeys the 
inelastic Boltzmann equation 



(at+vV)/ = J[/,/]+ij.p/, 



(1) 



where J\f, /] represents the collisional integral given 
by[nilI3 

J[f,h] = a"-^ J d& I dvie(g-^)g-<T 

x(e-Vi-I)/(r,v,t)/^(r,vi,t) (2) 

where Q{x) is Heviside's function, g = v — vi, and & as 
the unit vector along the line connecting centers of mass 
of contacting particles. The operator is the inverse 
of the collisional operator b which are defined as 



6g = g - (I + e)(g • 

1 + e, 



b g = g- 



(3) 
(4) 



where e is the coefficient of restitution which is ranged 
< e < I. Here we assume that e is a constant 
for the simplification of our argument, though the ac- 
tual coefficient of restitution depends on the impact 
velocity,!^ IT^ [l5l |. The effects of the impact veloc- 
ity dependence of e to macroscopic hydrodynamics can 
be seen in ref.^^. It should be noted that the oblique 
impacts have important contributions in actual inelastic 
collisions. ^3 El We, however, assume that the effects of 
inelastic oblique collision can be neglected, which may be 
justified when particles are smooth hard-core particles. 

The Fokker-Planck operator Lpp in eq.([3) represents 
the driven force coming from the Langevin force as 



L 



FP 



70 



d_ 
dv 



[V. 



(5) 



where the first term and the second term represent 
the frictional force and the thermal activation, respec- 
tively. In general, the temperature of the heat bath 
Tb is different from granular temperature T. If there 
is no contribution from the collisional integral, the dis- 
tribution function is relaxed to an equilibrium state as 



2 



f ^ feq oc exp[—mV'^/2TB]- The Fokker-Planck opera- 
tor may be categorized into one of the white noise 
thermostat L^hite — 7o7's/™(9^/c'v^). However, the 
contribution of Lpp to the distribution function, the 
transport coefficients and stabiHty of the homogeneous 
state is different from the pure white noise thermostat. 
We also note that the frictional force in the first term 
of (0 has the inverse sign of the Gaussian thermostat 7] 
Lcauss = / d"vyV which is also used in simulation 

of molecular hydrodynamics. ^3 

Hydrodynamic variables to characterize the macro- 
scopic behavior of the gas are the number density, the 
velocity field and the granular temperature defined by 



n(r,i) = j dvf{r,v,t), 
i(r,f)u(r,t) = y"dvv/(r,v,i). 



(6) 
(7) 



-n{v,t)T{v,t) = j dv-mV2/(r,v,i), (8) 

where V = v — u. The integral of the coUisional invari- 
ance multiplied by J[/, /] over V is zero. Since the loss 
of kinetic energy in each collision is given by 



1 - p2 



-r7i(g • 0-) 



(9) 



the following relation holds 



dvimt/V(/,/) = -^TC[/,/] (10) 



approximately 



)ling 



Here, the co oling rate C in ea. (|10|l can be evaluated 

C^C^") = i^^o(d + 2)(l + Aa2), (11) 

where = TT-^/'^na'^-^{T/m)^/'^md/{d + 2), 03 is the 
forth cumulant defined by 



a2 



d < > 
d + 2 < y2 >2 



- 1, 



(12) 
(13) 



which will be determined later. 

The balance equations for hydrodynamic variables are 



Dtn + nV • u = 0, 
DtUi + {mn)^^V -jPi-j 0, 
2 

dn 

= 27o(Ts-r), 



AT + — (i^yV,u, + V-q)+TC 



(14) 
(15) 

(16) 



where Dt = 9t -I- u • V. The pressure tensor and the 
heat flux q are respectively defined by 



(17) 
(18) 



P,, ^ mj dvy,VS/(r,v,i), 
q = dvy2v/(r,v,i). 



We adopt the Chapman- Enskog method'2'l, where 
space and time dependences appear through hydrody- 
namic variables. The expansion parameter is regarded 
as the magnitude of spatial inhomogeneity. Thus, we ex- 
pand / around the homogeneous solution /^"^^ as / = 
/^"^ +e/^^^ -l-e^/^^^ + • • • , derivative is also expanded as 
dt — dt^^"^ + edt''^^ -I- • • • . Here we apply the Chapman- 
Enskog method for dilute granular gases developed by 
Brey et al.^Bi l2lj| and SantosQ to driven systems. To 
remove the ambiguity of the distribution function, we 
impose the solubility conditions in which hydrodynamic 
variables are unchanged from the evaluation by (v, t). 



III. HOMOGENEOUS STATES 

As the first step of Chapman-Enskog method, we need 
also to obtain the homogeneous solution of the inelastic 
Boltzmann equation. We usually assume the scaling form 



/(v, t) - nva{t)-^f{c, r), c = \/va{t) 



(19) 



with dr — LOEdt and VQ{t) — ^J2T /m. The forth cumu- 
lant introduced in (|12|l for the scaling function is related 
to < c-* >= / dcc^fic) as < >= d{d + 2){a2 + l)/4. 
Here lue is Enskog's collision frequency given by 



d + 2 



4 V TT 

For the calculation we need to obtain 
^fe = - / dcc''J[fJ] 



(20) 



(21) 



The cumulants and fj,k can be evaluated by an approxi- 
mate expansion of Sonine polynomials^^- For example, 
/i2 is evaulated as 



^^"2(^-^)71^(^+16"^)' 



II4 is estimated as [ij 



A*4 



^ndiAi+a2"A2), 



(22) 



(23) 



with 



A2 



l-e^ 3 2^ 

^{l-e'){10d + 39 + 10e') 
+ i(l + e)(d-l). 



(24) 



It should be noted that these evaluations are based on 
two approximations: (i) the truncation of the first Sonine 
expansion, and (ii) the linearization of a2. If we adopt 
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the first assumption, the linearization of 02 gives a nice 
evaluation. |22| (iii) the direct comparison of the trans- 
port coefficients obtained by the Monte Carlo simulation 
and the linearized approximation gives good agreement 
in free cooling states. [21j However, nobody knows how to 
converge the Sonine expansion for driven granular gases. 
Thus, we may need to check the convergence of the So- 
nine expansion as in the case of elastic particles. |23j| 

Now, let us discuss the time evolution of temperature 
field. Equation 1)16(1 becomes 

dt^°^9 = 270 - (270 + CW (25) 

in the homogeneous state with the Langevin thermostat, 
where 6 = T/Tb- Assuming = 9^^^ + a29^^^ + 0{a2^), 
ea. (P5|l is reduced to 

5,0(0) =2{7-(7 + C)^(")} 



(26) 

in the lowest order, where 7 = Jo/uje and ( — (1 — e^)/2d. 
This equation has the solution 



QiO) 



(e(°)(o) 



-2(7-|-C)t 



(27) 



with 9oo = 7/(7 + 0- 

In the scaling limit, the inelastic Boltzmann equation 
(P) is reduced to 

%a./,c,.,.i[/,/K(l-^)|.,c;,c,.), 



7 ~ 



(28) 



where J[f,f] = n'^{a/vo{t))'^^^J[f,f]. From the equa- 
tion for < c'^ >= d{d + 2)(a2 + l)/4 we obtain the equa- 
tion of lowest order of 02 : 



dra2 



4C-ii 



.19; 



4j-i2] + Jc«2'(29) 



1^ 

0(0) 



A. 



(30) 



with Ai = 8Ai/d{d + 2) and A2 = 8A2/d{d + 2). The 
solution of H3UI) is given by 



02 (r) 



(02(0) -a2' 



a2 

xexp[(-C-A2)r-47y^ ^ 



where 



02^ 



Ai-4C 



fC-4(7 + C)-A2 



(31) 



(32) 



Although we present the result of linearization of 02 , it is 
possible to obtain the exact steady values of 9 and 02 for 
(|25|l and H29|l. The result is presented in Appendix, but 
the differences between the exact values and the result 
from linearized approximation are invisible (Fig.l). We 
also compare 02 in ea. H31|) and 9 in (|27|l with the numer- 
ical solution of H29|) and (|25|l to confirm the validity of 
the linearization of 02. 
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FIG. 1: Steady 02 as a function of e for d = 3 and 7 = 0.1. 
The solid line and the solid circles represent the result of 
linearized approximation and the exact one, respectively. 



IV. THE DETERMINATION OF THE 
TRANSPORT COEFFICIENTS 

Here, we explicitly obtain the result of the transport 
coefficients for heated granular gases. The solution /(") 
is isotropic so that the zeroth order pressure and the heat 
flux are given by 



q(0)=0 



(33) 



where p = nT is the hydrostatic pressure. 

The flrst order equation of the Boltzmann equation 
becomes 

(5,(")+L-L^p)/(i) = -(a,W+v.V)/(") 

= -(A^'^ + V- V)/(") (34) 

with Dt'^'^ = dS^^ -I- u • V. Here the linear operator L in 
(EU is defined by 

i/(i)^_j[/(o)^j(i)]„j[j(i)^j(o)]_ (35) 

Multiplying both side of H34|) by mViVj and integrate 
over V, we obtain 



(9t(°) -I- iy)P,/'^ + n,/i) = -pX.kiVkUi (36) 



where 



H,/i) = -m 



J dVV,VjLFpf(-^\ (37) 



Aijfei = SikSji + SiiSjk — —SijSki- (38) 
The solution of H36() can be written as 

PyW = -r;A,,fc,VfcM, (39) 
where ry is the viscosity. 



4 



It is possible to obtain v in ea. (|36|) through the relation 



(40) 



The evaluation of v is independent of the existence of the 
thermostat. We have evaluated v as IB: 



ly 3 



<^--4;^(l-e + ^d)(l + e)(l-la.). (41) 
With the aid of 



dt^''P..'''-[^-lo-^^)p./'\ (42) 

eqs.lO,® and ^ lead to 

7] 27* + 1 



V 



770 7*(l + ^"') + '^r,*-C*/2 



(43) 



where 7* = 7/1^0, i^tj* — v jv^) and C* = C/^'o- %(7) is the 
viscosity for e = 1 which is different from the value r\e of 
the elastic gas as — 1^0(2^* + 1). The steady value of 
77* is obtained when we substitute 0°° and a2°° into gSJ. 

Let us consider the heat flux. Multiplying both side of 
ea. (|34|l by mV'^'V/2 and integrate over V we obtain 

(a/°)+z.')q'i)+Q« = _^(i + 2a2)^VT 

2 m 

-^«2-Vn, (44) 
2 m 



where 



(45) 



Here, u' has already been calculated as|21 
, u' l + e,d-l 3 



:[^ + -(d + 8)(l-e) 



d ' 2 16 
4 + 5rf-3(4 - d)e 
512" 



-a2. 



The heat flux is described by 



(46) 



(47) 



where 77 and k are the shear viscosity and the thermal 
conductivity, respectively. The other transport coeffi- 
cient /i appears only is granular gases. Through the sub- 
stitution of and lUni) with the result of Q^^^ into iB^ 
we can obtain k, and /i. 

From the scaling form H19|l the following relations 
should be satisfied: 

9t(")q(i) = (2C + 37o-7o/f?)^VT 

+{m[^C + 370(1 - e-')] -t- ^}Vn(48) 
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FIG. 2: The time evolution of rf (solid line), tC (dashed 
line) and /i* (dot-dashed line) for e = 0.9, d = 3, 7 = 0.1 and 
61(0) = 2. 



On the other hand, 11^^^ and Q*^^^ have the relations 



(49) 



Thus, K and can be obtained from (^511,63) and ll^ 

as 



d + 2 H- 2a2 nT 
2 1^' -2C + jo/0~m' 

n ^ 2 "2 „ 



ly'-U-^lo/O' 



(50) 
(51) 



Thus, K* — k/kq which is the normalized heat conduc- 
tivity by its value of e = 1, ^0(7), is given by 



l + 2a2 



, d~l 

^ d ''i^x* -2c* + r/d 

K* a2{vn*+l*/e) 
1 + 2a2 Vn* - |C* + 37*/^ 



(52) 
(53) 



Their steady values are evaluated replacing 9 and 02 by 
their steady values. 

Figure 2 shows the time evolution of 77* , k* and fi* for 
e = 0.9,c? = 3,7 = 0.1. Figure 3 shows the steady values 
of 77*, K* and /i* as functions of e for d = 3 and 7 = 
0.1. As we can see from figures, /i* can be comparable 
with others. This situation has not been realized for free 
cooling systems. 

It is also remarkable that /i* is not proportional to 
02 and keeps positive for the region of 02 < even in 
the steady limit. This result is contradicted with the 
framework by SantosQ. The difference comes from the 
following. In our system, the system is heated by Tb uni- 
formly, but the temperature field T contains spatial fluc- 
tuations. Thus, dt^'^'^c^^^ cannot be zero even in steady 
states, because this contains the term in proportion to 
V[27oTb + (27o + C)T]. 
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reported that deviation of / from the Gaussian in driven 
granular gases is large. So we may need to check the 
convergence of the Sonine expansion and to compare the 
theoretical prediction with simulations. We also need to 
look for the possibility to apply our result to explain the 
data in actual experiments for suspensions. 

The author would like to thank A. Santos for fruit- 
ful discussion. This study is partially supported by the 
Inamori Foundation. 



APPENDIX A: STEADY VALUES OF aa AND 6 



FIG. 3: The steady values of rj* (solid line), k* (dashed line) 
and /i* (dot-dashed line) for e = 0.9, d = 3, 7 = 0.1. 



V. CONCLUDING REMARKS 

We have derived hydrodynamic equations based on a 
systematic Chapman-Enskog method for dilute granular 
gases driven by the Langevin thermostat in this paper. 
We have determined all of the transport coefficients rj, k 
and fi appear in Navier-Stokes order as a function of the 
restitution coefficient e. 

The result is based on the linearized approximation 
of 02 in the first order truncation of Sonine expansion. 
Although we believe that this approximation gives a nice 
evaluation, nobody knows its theoretical background and 
quantitative validity. In particular, Pagnani et aL^lOj has 



As mentioned in the text, it is possible to obtain exact 
steady solutions of (|25() and H29(l . The result is so com- 
plicated and difference between them and the linearized 
solutions is small that we do not use the exact form for 
later discussion. Here we present the exact solutions: 



a2- 



16(1 - 3e2 + 2e^ 



(Al) 



where Qi = 73 - 32e - 75e^ - 30e'^ + 64^^ + 8d(7 + 4e - 
3e^ + I67), and 

f) ^ .i2^ 

°° Q2 + 64d^^^ + Q3-dil + e)Qi' ^ ' 

where Q2 = 2(1 + ef{19 - 46e + 336^ - 126^ + 6e^), 
Q3 = 8^27(11 + 4e - 7e2 -f- 167) and Q3 = -28 + 6^(28 - 
3O7) - 1377 + 6e3(57 - 2) -I- e(12 + 169j). 
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